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ABSTRACT 

We present a clear detection of a broad Lya absorber (BLA) with a matching Ovi hne in the 
nearby universe. The BLA is detected at z(Lya) = 0.01028 in the high S/N spectrum of Mrk 290 
obtained using the Cosmic Origins Spectrograph. The Lya absorption has two components, with 
6(11 i) = 55 ± 1 km s~^ and 6(H i) = 33 ± 1 km s^^, separated in velocity by w ^ 115 km s^^. The 
Ovi, detected by FUSE at z(Ovi) = 0.01027, has a &(Ovi) = 29 ± 3 km s'^ and is kinematically 
well aligned with the broader H i component. The non-detection of other ions such as C ii. Si ii, 
Fe II, C III, Si III, C IV, Si iv and Nv at the same velocity as the BLA and the Ovi implies that 
the absorber is tracing highly ionized gas. The different line widths of the BLA and Ovi suggest a 
temperature of T = 1.4 x 10^ K in the absorber. Photoionziation, collisional ionization equilibrium as 
well as nonequilibrium collisional ionization models do not explain the ion ratios at this temperature. 
The observed line strength ratios and line widths favor an ionization scenario in which both ion- 
electron collisions and UV photons contribute to the ionization in the gas. Such a model requires a 
low-metallicity of ^ —1.7 dex, ionization parameter of log U ~ —1-4, a large total hydrogen column 
density of 7V(H) - 4 x 10^^ cm'^, and a path length of - 400 kpc. The line of sight to Mrk 290 
intercepts at the redshift of the absorber, a megaparsec scale filamentary structure extending over 
^ 20 deg in the sky, with several luminous galaxies distributed within ~ 1.5 h^^ Mpc projected 
distance from the absorber. The coUisionally ionized gas phase of this absorber is most likely tracing 
a shock-heated gaseous structure, consistent with a few different scenarios for the origin including 
an over-dense region of the WHIM in the galaxy filament or highly ionized gas in the extended halo 
of one of the galaxies in the filament. In general, BLAs with metals provide an efficient means to 
study T ~ 10*^ — 10^ K gas in galaxy halos and in the intergalactic medium. A substantial fraction of 
the baryons missing from the present universe is predicted to be in such environments in the form of 
highly ionized plasma. 



1. INTRODUCTION 

Cosmic baryon census estimates have shown that in 
our present universe only ~ 6 — 10 % of the primor- 
dial baryons are entrapped in galaxies. The bulk of 
the baryons have not collapsed into luminous structures 
but are in circumgalactic regions, and in the unvirial- 
ized l arge scale inte rgalact ic filaments (iM ulchae yet all 
[l99l iFukugita et a l. 1998; Fukugita fc Pe ebles 2003). A 
significant fraction (~ 50%) of these baryons outside of 
galaxies has already been discovered by STIS / HST and 
FUSE UV spectroscopic observations of the z < 0.5 
intergalactic medium (IGM). The low-z IGM exists as 
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a mix of multiple gas phases. The photoionized diffuse 
IGM with T < 10** K traced by the Ly a forest accounts 
for 30% o f the baryonic mass density (Shull et al. 199 
Penton et al. 2004; Lchner et al. 2006; Danforth & Sliul 
2008fl . The population of highly ionized metal line 



absorbers, notably Ovi systems discovered closer to 
galaxies, possibly contain an additional ^ 10% of the 
mass-fraction if their ionization is dominat ed by colli- 
sions o f electrons with ions in the plasma (iTripp et alj 
'2000bi; 'Richter et al.' '200 l IDanforth fc Shulll 120051: 
j Savage et al.n2005i: ,Stocke et all 120061: IDanforth fc Shulll 



120081: ITHto et al.H2008HNaravanan et al.ll2010[ ). The re- 

maining ~ 60% of the baryonic mass-fraction in our 
present universe remains to be discovered. Hydrody- 
namic simulations of structure formation predict that 
this missing portion from the baryon inventory is in 
a gravitationally shock-heated phase with temperatures 
in the range of T '-^ 10^ — 10 ^ K and densities of 
nn (0.1 - 10) X 10-5 cm-3 dCenfc Ostriked [1991 
IDave et al.ll2001l: IValageas et al.ll2002D . The simulations 
convey that the reservoir of cooler intergalactic gas gets 
heated to warm - hot temperatures as a result of shocks 
produced during the gravitationally collapse of matter 
to form galaxies and clusters in the over-dense regions of 
the cosmic web. Discovering this warm-hot intergalactic 
medium (WHIM) is one of the main science drivers for 
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the C osmic Origins Spectrograph (COS) on HST (|Shulll 
[2009h . 

Currently, UV absorption line spectroscopy is the most 
promising observational technique for detecting the com- 
ponent of the WHIM with T - 10^ - 10^ K. Among 
metal lines, Ovi AA1031,1037 A are the most sensi- 
tive for observations of high ionization gas, due to the 
high cosmic abundance of oxygen and the large values 
of oscillator strengths of the 1032 A and 1038 A reso- 
nant lines. Under collisional ionization equilibrium con- 
ditions, Ovi ionization fraction reaches its peak value 
of /(^^^^ ^ iV(Ovi)/iV(0 ) = 0.187 at T = 3.2 x 10^ K 
(IGnat fc Sternberg! [2nn7h . Thus, Ovi can be a tracer of 
the warm phase (i. e. the lower temperature phase) of the 
shock-heated gas (iTripp et al.l l200l1: iSavage et all 120021: 
iDanforth &: Shullll2008[ ). Theoretical estimations predict 
~ 20 — 30% of the missing baryonic matter is traced 
by Ovi absorpt ion systems with VKr(1032) > 20 mA 
()Cen et al.ll2001[ ). However, observationally confirming 
the value of rib(Ovi) has been a challenge. Detectable 
amounts of Ovi can also be produced in low-density 
gas (rijj < 10~^ cm~^) photoionized and heated by a 
hard radiation field, in which case the absorber may 
not be tracing a large-scale struc ture (ISavage et al.lll998l: 
Prochaska et all [200 4: Thom & ChenI 120081: iTripp et al.l 
20081 lOppenheimer fc Dava,2009ir The current sample 
of known Ovi systems seems to suggest that they are 
a heterogeneous mix, with evidence for photons domi- 
nating the ionization in some cases, and electron colli- 
sions in others. Detailed ionization modeling with tight 
constraints on column densities, Doppler widths and ve- 
locities of H I, Ovi and other metal ions are necessary 
to understand the ionization mechanism in each absorp- 
tion system and its physical origin (Savage et al. 200j; 
iNaravanan et al.ll2010D . 

An alternative to using highly ionized metals is to 
search for absorption from the trace amounts of atomic 
hydrogen (ionization fraction, /jj j = iV(H i)/A^(H) < 
10~^) present in the warm - hot gas. At T > 10^ K, 
an absorber with a total hydrogen column density of 
A^(H) ^ 10^^ cm~^ will produce Lya absorption with 
column densities of A^(H i) < 10^'' cm~^ and Doppler 
widths of &(H i) > 40 km s^^. Accurate measure- 
ments of such broad lines in the quasar spectrum re- 
quire high S/N , well-determined continua, and a good 
understanding of the fixed-pattern noise. Broad Lya 
absorbers (BLAs) in the z < 0.5 universe were discov- 
ered in the existing ar c hive of HST /STIS quasar spectra 
(iSembach et al.l 12004 iRichter et all 120061: iLehner et ahl 
120071 IDanforth et al.ll2O10ir The most serious systematic 
uncertainty associated with these BLA measurements is 
whether the observed broad line width of Lya is domi- 
nantly thermal. Due to the shallow and broad shape of 
the absorption profile, particularly when the H i column 
is low, it is often difficult to rule out, from moderate S/N 
spectra, the possibility of line blending or turbulence con- 
tributing to line width. High S/N spectroscopy is there- 
fore crucial to confirm the thermally broadened Gaussian 
absorption profile characteristic of BLAs. At low-z BLAs 
account for a baryonic mass-densit y fraction sim i lar to 
galaxies or possibly even greater ([Richter et al.l 120061 : 
IDanforth et al.ir2010|) . Thus discovering more BLAs is 



central to the continuing search for the reservoirs of miss- 
ing baryons in the present universe. Observations of the 
low-z universe also opens up the possibility to search for 
an answer to the allied issue of how these baryons are dis- 
tributed in relation to galaxies in large-scale structures. 

In this paper, we report on a BLA in the very nearby 
universe detected in the high S/N spectra of Mrk 290 ob- 
tained with the HST/ Cosmic Origins Spectrograph. The 
BLA is at the same velocity as an O vi absorber identified 
in the FUSE spectrum of the same sight-line. The BLA - 
O VI system is tracing collisionally ionized gas associated 
with a large-scale galaxy filament. We describe in detail 
the COS observations (sec 2.1), the properties of the ab- 
sorber (sec 3), and the details on the galaxies in the fore- 
ground field (sec 4). The physical properties of the ab- 
sorber are investigated using the predictions made by dif- 
ferent ionization scenarios (sec 5). Based on the observed 
properties and the ionization analysis, we speculate on 
the origin of the absorption system and its baryonic con- 
tent (sec 6). The cosmology that we adopt is ACDM 
with = 0.3, r^A = 0.7 and a dimensionlcss Hubble 
parameter of /i iJo/(100 km s^^ Mpc"^) = 0.71. 

2. OBSERVATIONS 

2.1. COS Observations 

The design capabilitie s of H ST / CO S are described in 
detail by Green (20011): iFroning fc G reen (200 l) and in 
the u pdated COS Instrument Handbook (Dixon et al.l 
20101). The infii g ht pe rfor mance of COS is discussed by 
Osterman et al.l (|2010l ) and in the numerous instrument 
science reports on the STScI COS websitcQ.The COS ob- 
servations of Mrk 290 were carried out on 28 October 
2009 as part of cycle 17 COS GTO program ID 11541 
(P.I. James Green). The details of the observations are 
listed in Table 1. The G130M and G160M COS gratings 
were chosen to obtain spectra in the far ultraviolet bands. 
The G130M observations were split into four exposures of 
^ 964 second integrations at grating central wavelength 
positions of 1291 A, 1300 A, 1309 A and 1318 A. Sim- 
ilarly, the G160M observations were composed of four 
~ 1200 second integrations with central wavelength set- 
tings of 1589 A, 1600 A, 1611 A, and 1623 A. Ah expo- 
sures were obtained in the time-tagged, photon-counting 
mode. The different grating central wavelength configu- 
rations resulted in photons of the same wavelength get- 
ting recorded at different regions of the detector which 
subsequently reduced the effect of fixed pattern noise in 
the final coadded spectrum. This set up also helped to 
cover the 18 A wavelength gap between segments A 
and B of the FUV detector for the various central wave- 
length settings. 

The reduction of the raw frames into ID spectra 
was carried out using the COS pipeline software Calcos 
(v2.11) with the appropriate calibration reference files. 
The one-dimensional science data files generated by the 
pipeline had a dispersion scale of 0.01 A corresponding to 
6 wavelength bins per resolution element. The individual 
integrations were coadded weighted by their respective 

^ http:/ /www.stsci.cdu/hst/cos/documents/isrs 
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exposure times. The final spectrum had a wavelength 
coverage from 1136 A to 1796 A. The signal-to-noise ra- 
tios (per 18 km s~^ resolution element) of this combined 
final spectrum at 1250 A, 1450 A, and 1650 A are 35, 24, 
and 16 respectively. 

Ghavamian et al. (2009) have determined the resolv- 
ing power of the spectrograph by detailed modeling of 
the line spread function (LSF) at various wavelengths. 
The spectral resolution is found to be wavelength depen- 
dent with values in the range A/AA ~ 16,000 - 21,000 
for the G130M and G160M gratings, where AA refers to 
the width at half strength of the LSF which has broad 
wings containing 25% of the LSF area. The resolution 
is maximum at near-UV wavelengths and declines mono- 
tonically towards lower wavelengths. 

To determine zero-point offset in the COS spec- 
trum we compare the velocity centroids of the wlsr = 
-138 km s^i ISM absorption lines of Sii A 1254 A, 
Fe II A1608 A a nd Al ii A1671 A from Complex 
C (jWakkeii I2OO60 . These lines appear at wlsr — 
— 159, —154, and —153 km s~^ respectively in the COS 
spectrum. A velocity correction of -1-17 km s""'^ was 
therefore applied to the COS data, which is approxi- 
mately one resolution element. We estimate residual er- 
rors of ~ 5 km s~^ (1 cr) in the velocity calibration. 

2.2. FUSE Observations 

The properties of the FU SE satellite and its p erfor- 
mance are described in detail by IMoos et all (j200Q) and 
ISahnow et all (|2000D . The processing of the FUSE data 
for this sight line is described in great detail elsewhere 
(jWakker et all [20031: [Wakker 2006). We only provide a 
short overview here. The archived FUSE observations 
for Mrk 290 were obtained over the time period from 2000 
to 2007 (exposure IDs: P1072901, D0760101, D0760102, 
E0840101, E0840101, E0840102). The reduction of these 
data sets was carried out using version 2.4 of the FUSE 
calibration pipeline. The shifts in wavelength and cen- 
tral velocities in the individual reduced exposures were 
corrected by aligning the Milky Way ISM lines with the 
LSR velocity of the 21 cm H i absorption from Com- 
plex C (wLSR = —138 km s~^). Once the velocity shifts 
were applied to each observation, the data in the LiFlA 
and LiF2B channels of the detector were combined to- 
gether to produce the final spectrum. 

3. PROPERTIES OF THE ABSORBER 

At a heliocentric velocity of whelio — 3081 km s^^, 
the z(Ovi) = 0.01027 absorber is one of the nearest 
extragalactic Ovi absorption systems. It was discov- 
ered by Wakker & Savage (2009) in the FU SE spectrum 
of Mrk 290 through the identification of Ovi A1032 A 
line. The detection was claimed as uncertain since the 
red member of the Ovi doublet is confused with the 
Galactic An A1048 A line. A hint of the Ovi A1038 A 
came from the extra absorption wing seen at positive 
velocities of the Galactic An line. The Lya for this 
absorber was outside the FUSE spectral window, and 
the Ly/3 at this velocity suffers from Galactic Cn con- 
tamination. Nonetheless, 1 Wakker fc Savaga ((2009|) re- 



ported the > 3 cr absorption detected at A 1042.5 A as 
Ovi A1032 A after ruling out the possibility of the iden- 
tified line being interstellar absorption or higher order 
Lyman lines along this sight line. The detection of Lya 
in the C OS G130M spectrum at z (Lya) = 0.01028 con- 
firms the I Wakker fc Savaei (|2009l ) identification of this 
absorption system. 

In Figures la - lb we show the continuum normal- 
ized profiles of Lya and Ovi on a velocity scale cen- 
tered at z(Ovi) = 0.01027. Also displayed are regions 
of the spectrum where absorption from other prominent 
ionic species are expected. The Lya profile is spread 
over a velocity interval of Au = 400 km s~^ in the rest 
frame of the absorber. The high S/N COS data show 
two distinct components in H i, separated in velo c ity by 

V - 115 km s-^ Using the iFitzpatrick fc Spitzeil ()1997[ ) 
software, we fit Voigt profiles to the Lya feature. The 
fit convolves the model profile with the G130M instru- 
mental spread function at the observed wavelength of the 
Lya feature. It is important that the specific line spread 
function determined for COS be used to fit the profile in 
order to minimize the impac t of the non-gaussian w ings 
on observed line profiles fsee iGhavamian et al.ll20d9l for 
a detailed discussion). 

The fit that best reproduces the observed ab- 
sorption profile is a two component model with 
log [7V(H i) cm-2] ^ 14.35 ± o.Ol dex and 6(H i) = 
55 ± 1 km s^^ at u = 2 ± 5 km s~^ and 
log [Ar(H i) cni-2] ^ 14.07 ± 0.01 dex and 6(H i) = 
33 ± 1 km s^^ at w = 116 ± 5 km s^^. This model 
fit is displayed in Figure 2. The fit reveals some weak 
residual absorption in the positive and negative velocity 
wings of the Lya. A broader, shallower H i component 
could explain this excess absorption in the wings. How- 
ever, because of blending with the stronger H i compo- 
nent at w = 2 km s~^, the fitting routine is not able to 
yield meaningful parameters for an extra component at 

V ^ —125 km s^^. The contribution to the total column 
density from such a shallow component is going to be 
marginal {Wr 30 mA, log [7V(II i) cni"^] < 13.0) when 
compared with the contributions from the two stronger 
components at w 2 km s~^ and v ~ 116 km s~^ re- 
spectively. 

The 6 = 55 km s~^ Doppler width for the v = 2 km 
H I component suggests that this is a broad Lya ab- 
sorber (from here on, wc often refer to this component in 
H I as the BLA). If the line is purely thermally broad- 
ened, then 6 = 55 km s^^ line width implies a gas tem- 
perature of T = 1.8 X 10^ K. The most significant sys- 
tematic uncertainty that could affect a 5-measurement 
is blending of closely separated kinematic components. 
The residual absorption in the blue wing of the BLA 
suggests some kinematic complexity. However, the BLA 
and the v ~ 116 km s~^ components of the Lya profile 
appear highly symmetric suggesting very little infiuence 
of this low column density component on the measured 
b values. 

We take note of the fact that the Lya profile could 
be subjected to saturation that is unresolved by COS 
(FWHM ~ 17 km s~^). The presence of instrumen- 
tally unresolved saturation becomes evident when appar- 
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ent column density profiles of different electronic transi- 
tions from the same energ y level in an ion are compared 
(jSavage fc Sembacblll991[ ). Such an analysis is not feasi- 
ble in the case of the z = 0.01027 absorber as the Ly/3 
covered by FUSE is strongly affected by contamination 
from Galactic ISM C ii absorption. The Ly7, which is 
^ 18 times weaker compared to Lya in line strength 
(/A), falls at the edge (A - 982.5 A) of the FUSE cov- 
erage where the quality of the data is poor {S/N ^4 — 6 
per resolution element). Even so, the H i upper limit 
on column density of 10^'^ "'^ cm~^ given by the low S/N 
Ly7 data suggests that there is no significant saturation 
in either component of Lya. 

The COS spectrum offers coverage of several impor- 
tant low and high ions including C ii, Si ii, Fe ii, Si ill, 
C IV, Si IV and Nv (see Figures la - lb). The expected 
wavelengths for all these ions show non-detections at the 
3 cr significance level. We therefore quote in Table 2 
upper limits on the column density for these nondetec- 
tions. The lack of any low ionization absorption is indica- 
tive of high ionization conditions in the absorber. The 
C III A977 A line is covered by FUSE, but falls at the 
edge of the spectrum where the S/N is poor (see Figure 
la). The line is not detected at > 3 tr significance. 

A Voigt profile fit to the Ovi A1032 A line yields 
a single component at u = 1 ± 5 km s~^ with 
log [iV(Ovi) cm-2] = 13.80 ± 0.05 dex and 6(0vi) = 
29 ± 3 km s-i. We use a Gaussian of FWHM = 
25 km s^-^ in this fit model to account for the FUSE 
instrumental smearing. The profile fit is displayed 
in figure 2. Using these fit parameters, we synthe- 
size the Ovi A1038 A line and superimpose the pro- 
file on top of the data in figure 2. The compari- 
son shows the uncontaminated positive velocity segment 
of the Ovi A1038 A absorption consistent with the 
Ovi A1032 A profile. We a lso use the apparent opt i- 
cal depth (AOD) technique of lSembach fc Savagd (|1992f ). 
to derive log [7Va(0vi) cm'^] = 13.82 ± 0.08 dex and 
6„(0vi) = 35 ± 8 km s'^ for the Ovi A1032 A line. Our 
AOD measurements are similar to the Wakker & Sav- 
age (2009) measurement of the same feature (see Table 
2). The ba is the second moment of the apparent col- 
umn density distribution, which when corrected for the 
FUSE instrumental spread gives a 6-value that is con- 
sistent with the profile fit measurement. 

In Figure 3, we plot the apparent column density 
profiles of Lya and Ovi. It is evident that the Ovi 
absorption is kinematically well aligned with the BLA 
(i.e., the w = 2 km s^^ H i component). The formal 
velocity difference between the line centers is Aw = 
u(Ovi) — v{Lya) = — 1 ± 5 km s~^. If the two ions 
are tracing the same gas phase, then the Doppler widths 
of &(H i) = 55 ± 1 km s'^ and 6(0vi) = 29 ± 1 km s'^ 
implies a temperature of T 1.4 x 10^ K with a non- 
thermal line broadening component of 6nt = 26 km s~^0. 
The derived temperature suggests that the v = 2 km s~^ 
H I component is by definition still a BLA (thermal 

^ We solve for the temperature assuming that the observed b 
is a quadrature sum of the thermal and non-thermal b values, i.e., 

bobs = y^FT^. 



5 > 40 km s ^), even though ^ 50% of the observed 
line broadening is non-thermal. 

The wavelength region of the Ovi A 1032 corresponding 
to the narrower v ^ 116 km H i component is con- 
taminated by a feature that is unrelated to this system. 
The wavelength region corresponding to the Ovi A 1038 
line shows no Ovi absorption associated with this H i 
component. The 6(11 i) = 33 ± 1 km s^^ implies an up- 
per limit on temperature of T < 6.6 x 10* K, which is 
consistent with photoionization. No other low or high 
ionization metal lines are detected at the velocity of this 
H I component. In Table 2, we separately list rest- frame 
equivalent width and column density upper limits for the 
various ions corresponding to the z; ~ 116 km s^^ Hi 
absorption. 

4. ABSORBER - GALAXY ASSOCIATION 

A search through the NED database reveals several ex- 
tended sources within 1000 km s~^ of the BLA - Ovi ab- 
sorber in the foreground field of Mrk 290. The left panel 
of Figure 4a shows the projected distribution of these 
galaxies within a ~ 16 Mpc region. The distribu- 
tion pattern suggests that the sight line is intercepting 
at the velocity uhelio = 3081 km s~^ of the absorber 
an over-dense filamentary structure in the cosmic web, 
with the galaxies organized around the axes of the fil- 
ament. The megaparsec scale filamentary structure ex- 
tends 15° on the sky. In Figure 4a, we also plot the 
position of other AGN sight lines in the field. Among 
those, Mrk 876 lies 9.7° ('- 7 /i"^ Mpc) further along the 
filament and shows Lya, Ly/3 and Ovi absorption at 
t^HELio = 3480 km s~^ most like ly associated with the 
filament (| Wakker k Savagd [2009f) . The PG 1626 -I- 554 
sight line also probes this filament at a separation of 7°. 5 
(^ 5.6 Mpc) in the sky in the opposite direction to 
Mrk 876, but does not show any H i absorption near the 
expected velocity. The other AGNs shown in the Mrk 290 
field (e.g. Mrk 279, Mrk 817, PG 1351 + 640) lie farther 
away from the filament and do not have any metal line 
or H I absorption detected near whelio ^ 3100 km s~^. 

The Mrk 290 field is inside the SDSS footprint. The 
galaxy sample in this region of this sky is nearly 100% 
complete down to a magni tude of g = 17.5 which cor- 
responds to - 0.03L* (Bla nton et all 120031) . From the 
more than two dozen galaxies within a projected dis- 
tance of ~ 1.5h~^ Mpc from the absorber, NGC 5981, 
NGC 5982, NGC 5985, NGC 5987 an d NGC 5989 are 
mem bers of the loose group GH 158 (IGeller fc Huchral 
Il983f ). The galaxy distribution shows at least nine other 
galaxies within ~ 1.5 Mpc impact parameter and 
At; = 300 km systemic velocity with respect to the 
absorber, but smaller in size and brightness compared to 
the GH 158 group members. Some of the other brighter, 
larger galaxies along this filament are also members of 
other galaxy groups (e.g. GH151). 

From all these known extended foreground sources, 
the luminous galaxy NGC 5987 is closest in impact 
parameter to the absorber. NGC 5987 [a(J2000) = 
15''39"57.4^(5(J2000) = +58°04'46"] is at a projected 
distance of p = 424 kpc from the sight line. 

The galaxy has a redshift oi z — 0.0100 (vhelio = 
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3010 km Ide Vaucouleurs et al.l 11991] ) which gives it 
a |Ai;| — 71 km systemic velocity with respect to 
the BLA - Ovi absorber. The galaxy has an estimated 
isophotal major diameter of D25 = 51.7 kpc, and is 
one of the la rgest spiral galaxies in the nearby universe 
(iRomanishirJ [l983i) . Although the impact parameter is 
large, the small velocity separation and the large size 
of the galaxy makes it possible for the absorber to be 
in some way related to NGC 5987's circumgalactic en- 
vironment. Such an inference is consistent with sta- 
tistical results from a number of past absorber-galaxy 
surveys where > 0.25L* gala xies were found within 
500 kpc of Oyi absorbers (IiHdd fc Savaed l2000al: 
Tripp et al.ll2000bt iTumlinson et all 120051: iStocke et alJ 
20061 ICoo ksev et al.' '2008': 'Chen fc MulchaevI l200a 
Lehner et a l. 2009; Wakkcr & Savage 200^ and also with 
results fro m cosmological simula tions of galaxies and Ovi 
absorbers (jGangulv et al.ll2008D . 

NGC 5987 is a spiral galaxy of Sb morphology with a 
B-band magnitude of = 1 2.72 or Mb -5 log h = -20.4 
(|de Vaucouleurs et al.lll99ll ). This corresponds to a lu- 
minosity of L ~ 2 L* assuming = —19.6. No pub- 
lished spectrum is available for this galaxy. In Figure 
4b, we display the orientation of the galaxy with respect 
to the quasar sight line. The large impact parameter 
of the absorber with respect to the galaxy disk results 
in low surface flux of UV photons from the galaxy in- 
cident on the absorber. Thus, we do not expect the 
galaxy's radiation field to significantly influen ce the ion- 
ization of hydroge n or metals in the absorber (|Fox et al.l 
|2005; Narayanan et"alll2010f ). NGC 5987 also does not 
appear in the ROSAT all -sky survey bright source cat- 
alogue ()Voges et al.l 119990 and therefore is unlikely to 
harbor an active AGN that could photoionize the sur- 
rounding intergalactic gas. 

Given the significant excess of galaxies with whelio 
3080 km s~^ in the field around Mrk 290, and the pen- 
cil beam size region probed by the BLA - Ovi system 
(whelio = 3081 km s~^), it will be difficult to conclude 
whether the absorption is tracing the gaseous halo of one 
or the other galaxy in the filament, or an intergalactic gas 
cloud that is part of the large-scale filament. We return 
to this discussion in Sec 5. 

5. IONIZATION CONDITIONS IN THE ABSORBER 

Understanding the ionization, physical conditions and 
metalicity in the absorber is crucial to discriminate be- 
tween the different possible explanations for its astro- 
physical origin. We explore a range of parameter space 
in both photoionization and collisional ionization mod- 
els and assess whether the models are physically realistic. 
We focus our modeling on the u ~ 2 km s~^ BLA compo- 
nent of the H I absorption which has O VI detected at the 
same velocity. No metals are detected for the other H i 
component at v ^ 116 km s"""^ and therefore useful con- 
straints on the physical parameters cannot be obtained 
for this component from ionization modeling. 

5.1. Photoionization 

The T = 1.4 X 10^ K temperature implied by the dif- 
ferent BLA and O Vi line widths is much higher than the 



kinetic temperature expected for a purely photoionized 
gas. Nonetheless, to elaborate on the difficulties asso- 
ciated with photoionization in explaining the origin of 
these lines, we investigate some best-fit photoionization 
models and their predictions. We use the ph otoioniza- 
tion code Cloudy [ver.C08.00, iFerland et all g,998)] to 
solve for equilibrium models that reproduce A^(Ovi) and 
the limits on all the other ions. The source of ioniza- 
tion is assumed to be dominated by the extragalactic 
background radiation at z = 0.01027 whose shape and 
intensity is as modeled by (jHaardt fc Madau 2001). The 
UV background spectrum used has contributions from 
quasars and from the population of high-z young star 
forming galaxies. The photoionization models are calcu- 
lated for different ionization parameters log OQ, metal- 
licity and the observed 7V(H i) = lO^^ -^^ cm~^ in the 
broad H i component a.t v — 1 km s~^. The iV(Ovi) = 
j^qIS.so q]j^-2 g^j^^j |.jjg lijjiiting column density ratios of 
N{C iv)/iV(Ovi), 7V(Nv)/Af(0vi), N{Si iv)/iV(Ovi) 
and A^(C iii)/A^(Ovi) serve as additional constraints for 
the models. In the models, we assume the revised so- 
lar elemental abundance ratios of [C/H]0 — -3.57 dex, 
[0/H]q = -3.31 dex, [N/H]0 = -4.17 dex and [Si/H]© = 
-4.49 dex given by Asplund et al. (2009) 

Shown in Figure 5 are the photoionization model pre- 
dictions for column densities of the various ionic species, 
as a function of the ionization parameter. The observed 
A^(C iv)/iV(Ovi) and iV(Nv)/A^(Ovi) limiting ratios re- 
strict the ionization parameter to log U > —0.6, corre- 
sponding to njj 6 X 10~^ cm~^ and metallicity to -1.8 
dex of solar. The A^(H i) = iQi^-SS ^^-2 jj 
dependent ionization fraction of /jj j < 10~^ implies an 
absorbing region with a total baryonic column density of 
A^(H) > 4 X 10^^ cm~^. The very low density and high 
total hydrogen column density results in a line of sight 
thickness of L > 3 Mpc for the absorbing structure. If 
the absorber is unvirialized, the large path length would 
result in lines with v > H{z) L ~ 210 km s~^ width due 
to Hubble flow, which is inconsistent with the observed 
H I and Ovi velocity widths. If the absorbing structure 
is well within the virialized halo of NGC 5987, it would 
have decoupled from cosmic expansion. The kinematic 
properties of the BLA would then largely depend only 
on the internal velocity field along the line of sight and 
the rotational velocity of the galaxy halo. It would be 
remarkable for a megaparsec scale absorbing region to 
not have any velocity field capable of inducing asymme- 
try or kinematic complexity in the line profiles. To the 
contrary, we find no such evidence in the observed line 
profiles of the BLA and Ovi which appear kinematically 
simple, symmetric and Gaussian. The predictions of the 
models are further complicated by the fact that the pho- 
toionization equilibrium temperature of ^ 0.6 x 10^ K is 
significantly smaller than the T = 1.4 x 10^ K given by 
the different BLA and Ovi line widths. Thus photoion- 
ization does not appear as a dominant process controlling 
the ionization state. We therefore turn to the alternate 
scenario of collisional ionization. 



^ Ionization parameter is defined as the ratio of the number 
density of photons with E > 13.6 eV to the hydrogen density, 
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5.2. Collisional Ionization 

5.2.1. Collision Ionization Equilibrium Models 

The temperature suggested by the observed hne widths 
is consistent with coUisional ionization conditions. First 
we consider the simplest scenario of colhsional ionization 
equihbrium (CIE) in which at a constant kinetic temper- 
ature the fraction of ions remain time-independent. In 
Figure 6, we compare observati ons with column density 
predictions from CIE models of ()Gnat fc Sternber32007l) 
for a metallicity of 1/10 solar. We find that CIE mod- 
els fail to explain all observed high ion column density 
limits simultaneously for acceptable temperatures. The 
predicted N{Oyi) is greater than iV(C iv) and iV(Nv) 
for T > 2 X 10^ K, which is higher than the temperature 
given by the line widths of the BLA and Ovi. Between 
the temperature range of T = 3.2 x 10^ K, where the 
Ovi ionization fraction peaks and T = 10^ K, the Ovi 
fraction declines rapidly by a factor of ^ 4 dex. Thus, 
at T = 1.4 X 10^ K, O.IZq metallicity and N{R i) = 
10^^-^^ cm"^ the CIE predicted iV(Ovi) is a factor 
of 3 dex smaller than the measured value and the 
iV(C iv) and iV(Nv) are significantly overproduced for 
solar relative elemental abundances. 

5.2.2. Hybrid of CIE & Ptiotoionization Models 

Even when conditions are favorable for collisional pro- 
cesses to dominate the ionization state, a realistic mod- 
eling of absorption systems should take into account the 
effect of UV metagalactic flux from quasar s and star 
forming galaxies on the absorber (iDanforth et al.l 120061: 
Richter et all 120061: iTripp et al.l m08i. iDanforth et al l 
20101 ). We therefore compute a series of hybrid models 
using Cloudy, allowing both collisional and photoioniza- 
tion reactions to contribute to ion fractions. The tem- 
perature in Cloudy is set to T = 1.4 x 10^ K given by the 
different BLA - Ovi line widths. The model that Cloudy 
would converge on will have the ionization and recombi- 
nation in equilibrium even though the heating and cool- 
ing rates are out of balance, due to the fixed temperature. 
In Figure 7 we show the column density predictions for 
these hybrid models at T = 1.4 x lO'^ K and different 
log U. The density of the absorber in these models is de- 
termined by the input number density of hydrogen ion- 
izing photons and the ionization parameter. The hybrid 
models that reproduce the measured iV(Ovi) and the 
column density ratios of Ovi with C iv and Nv require 

[Z/H] 1.9 dex, log U 0.6, njj - 4 x 10"^ cm"^ 

7V(II) ^ 1.8 X 10^° cm~^ and thus absorber path lengths 
L ^ 12 Mpc. The expected line width from Hubble 
broadening for such a large path length is incompatible 
with the observed H i and O vi line widths and therefore 
a CIE - photoionization hybrid model can be rejected. 

5.2.3. Non-equilibirum Collisional Ionization Models 

The temperature T ^ 10^ K implied by the data is 
in the regime where radiative cooling is efficient. The 
presence of metals can enhance the rate of cooling of gas 
such that significant departures of ion fractions from CIE 
will happen at these temperatures. The time dependent 
ionization fractions in a gas radiatively cooling from an 



equilibrium temperature of T 10^ K were calculated 
by Gnat & Stenberg (2007). In Figure 8, we illustrate the 
non-equilibrium collisional ionization (noneq-CI) curves 
for such a cooling gas flow using their model for metal- 
licites of Zq and O.IZq. In non-equilibrium cooling mod- 
els, the recombination process lags behind the declining 
temperature such that the ionization fraction for high 
ions will be higher compared to CIE for temperatures 
in the range 

104-5 _ 105-3 K. The effect is pronounced 
for higher gas metallicities. At solar metallicity and 
T = 1.4 X 10^ K, the Ovi ion fraction in noneq-CI mod- 
els is /(Ovi) = 2.63 x 10"^^ a factor of 1.4 dex larger 
than the Ovi ion fraction at CIE. From Figure 8, it can 
be inferred that for Z 0.3 Zq, the A^(Ovi) given by 
the model will be in agreement with the observed value 
for r = 1.4 X 10^ K. However, for the same model, the 
predicted N{G iv) ^ iV(Ovi) whereas observations sug- 
gest 7V(Ovi)/7V(C iv) > 8. The C iv prediction from 
this model would be consistent with the data only if the 
carbon to oxygen abundance ratio is reduced from solar 
by 1 dex. Thus, a simple noneq-CI scenario may not 
explain in a satisfactory manner the ionization in this 
absorber. However by including in these noneq-CI mod- 
els the influence of UV photons from the extragalactic 
background radiation, an acceptable ionization solution 
can be obtained. This is elaborated in the following sub- 
section. 

5.2.4. Hybrid of noneq-CI & Photoionization Models 

In Figure 9, we display the changes to Ovi and C iv 
fractions when collisions of ions with both electrons and 
photons are simultaneously influencing the ionization in 
the gas. The models in Figure 9 were computed using 
Cloudy for a fixed temperature of T = 1.4 x 10^ K given 
by the BLA - Ovi line wid ths and a fixed intensity of 
the lHaardt fc Madaul (|200lD (QSOs and galaxies) extra- 
galactic radiation field determined for the redshift of the 
absorber. The models are therefore a hybrid of equilib- 
rium photoionization and equilibrium collisional ioniza- 
tion reactions. The changes perceived for the ion ratios 
in these hybrid models can be used to estimate how the 
ion ratios would be affected in a noneq-CI cooling gas sce- 
nario when photoionization is included. It should not be 
forgotten that several broad assumptions are inherent in 
these ionization models (plane-parallel slabs of uniformly 
cooling gas with no internal temperature gradient). Few 
real absorbers would resemble closely with these condi- 
tions. We should only expect these computational mod- 
els to converge on a general acceptable solution rather 
that an exact fit to the full range of observed ionic col- 
umn densities. Our correction allowing for the effects 
of photoionization on the noneq-CI cooling gas is only 
approximate. A proper calculation would need to fol- 
low the changing ionizing conditions in gas cooling from 
T - 10^ K to T ~ 1.4 X 10^ K when simuhaneously 
subjected to the extragalactic radiation field. 

From the hybrid models shown in the left panel of Fig- 
ure 9, we draw the following conclusions. For low values 
of log U (i.e. high density), ionization by photons does 
not contribute significantly to the Ovi to C iv ratio. For 
log U < -4, iV(0vi)/7V(C iv) is similar to its CIE value 
ofiV(Ovi)/7V(C iv) 3.4dexatr= 1.4x10^ K. With 
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an increase in ionization parameter (i.e., a large number 
of ionizing photons per electron), photoionization begins 
to increasingly influence the ionization fraction of Ovi 
such that over the interval —4.5 < log U < —1.5, the 
iV(Ovi)/iV(H i) and iV(0vi)/7V(C iv) ratios are raised 
by a factor of ^ 3 dex. This region of the parameter space 
corresponds to the interval that separates models where 
Ovi is predominantly coUisionaly ionized from models 
where photoionization begins to predominate. Within 
this log U interval the C iv ionization fraction remains 
the same as CIE and begins to decline at log U > —2.5. 
We find that the change to the A^(Ovi)/A^(C iv) ion 
column density ratios seen in Figure 9 {left panel) is not 
sensitive to metallicity. This increase in the ionization 
fraction of Ovi from pure collisional values, due to the 
introduction of the ionizing radiation field can be used 
to interpret hybrid models of noneq-CI and photoioniza- 
tion. 

For [Z/H] = -2.0, the noneq-CI model predictions at 
T = 1.4 X 10^ K are iV(Ovi) - 10^°-^ cm'^, N{C iv) - 
10^^ cm^^, A^(Nv) ~ 10^^-^ cm~^ (refer right panel Fig- 
ure 9). The Nv prediction is already consistent with 
the upper limit given by the data, whereas C iv is over- 
produced by a small amount (0.1 dex higher), and Ovi 
underproduced by ^ 3.2 dex. In the hybrid model shown 
in the left panel of Figure 9, we can see that the O vi ion- 
ization fraction is enhanced from noneq-CI by ^ 3.2 dex 
for log U ^ —1, whereas the C IV ionization fraction is 
diminished by ~ 0.5 dex, making both column density 
predictions consistent with the measurements. However, 
the density njj ~ 10"*^ cm~'^, and total hydrogen column 
density of A'^(H) ^ 6 x 10^^ cm^^ corresponding to this 
ionization parameter implies a line of sight of ~ 2 Mpc 
through the absorbing region. Such a large path length is 
inconsistent with the observed Lya and Ovi line widths. 
Lowering the ionization parameter would decrease the 
amount of Ovi and increase the amount of C IV making 
the predictions disagree with measurements. 

For the hybrid model to agree with measurements, the 
log U has to be reduced (thereby narrowing the path 
length) simultaneously raising the metallicity in the col- 
lisionally ionized phase. For [Z/H] = -1.6, the iV(Ovi) 
will increase by ^ 0.7 dex from its [Z/H] = -2.0 value (see 
Figure 9) whereas the A^(C iv) will increase by ^ 0.4 dex. 
In the presence of ionizing photons the predicted 7V(0vi) 
will increase by an additional 2.5 dex making it equal 
to the measured value for log U ^ —1.5. For this hy- 
brid model, the A^(C iv) will be ~ 0.4 dex smaller than 
the noneq-CI predicted value of A^(C iv) ~ lO^'^-^ cm~^ 
making it approximately consistent with the 3 a up- 
per limit of 10^^'^ cm^^. The log U ^ —1.5 corre- 
sponds to rijj ^ 4 X 10^^ cm^^ and a total hydrogen 
column density of A^(H) ~ 4 x 10^^ cm~^ implying a 
path length of ~ 320 kpc for the absorbing region. Such 
a path lengt h is consistent with the L yg and Ovi line 
widths. The iGnat fc Sternberg! (pOOl noneq-CI models 
(shown in Figures 8 and 9) are available only for values of 
metallicity that differ by 1 dex. We have therefore used 
approximate scaling estimates for the noneq-CI column 
densities at [Z/H] = -1.6. 

There is a narrow range in metallicity and ionization 



parameter space where the hybrid model predictions will 
be consistent with the measured column densities and 
line widths. The values of [Z/H] = -1.6 and log U ~ —1.5 
are an upper limit on the metallicity and a lower limit 
on the ionization parameter, respectively. Increasing 
the metallicity would result in C IV being overproduced 
from the Ovi phase. For a lower metallicity of [Z/H] 
= -1.8, the log U ^ —1.3 in order to recover the ob- 
served 7V(0vi). This hybrid model would also be con- 
sistent with the A^(C iv) upper limit. However, the sug- 
gested path length is ^ 500 kpc. Thus [Z/H] = -1.8 and 
log U ^ —1.3 are therefore lower and upper limits re- 
spectively. Thus, within the framework of this simple 
noneq-CI and photoionization hybrid model, the BLA - 
Ovi gas phase has a metallicity of [Z/H] = — 1.7 ± 0.1 dex 
and ionization parameter log U = —1.4 ± 0.1 dex. 
The implied total hydrogen column density is very large 
(A^(H) ^ 4 x 10^^ cm~^) and therefore interesting. It is 
intriguing that even though the BLA is tracing collision- 
ally ionized hydrogen at T 1.4 x 10^ K, in the hybrid 
model much of the Ovi is produced through photoion- 
ization. 

It is important to emphasize here that the estimated 
metallicity and ionization parameter are dependent on 
the simplistic assumptions intrinsic to the ionization 
models and hence should not be considered as very ro- 
bust. For example, the predicted ^ 2% solar metallicity 
and njj 4 x 10~^ cm~'^ density would imply a radiative 
cooling time-scale of ^ 25 Gyr for the gas to cool from an 
initial temperature of T '-^ 10^ K 0. However, this time- 
scale estimation depends on the density and metallicity 
of the gas, and the initial temperature from which the 
gas begins cooling, all of which are not well determined 
for this absorber. A lower ionization parameter (higher 
density), higher metallicity or lower initial temperature 
would all result in lower cooling time-scales values that 
are consistent with the Hubble time-scale. 

The above hybrid model also has the limitation that 
it does not take into account the enhanced metal line 
cooling due to the increase in line emissions from excess 
Ovi produced by photoionization. A ~ 3 dex increase 
in the Ovi column density could lower the cooling time 
scales thereby altering the ionization states of C IV and 
Ovi. a more accurate treatment of the hybrid model 
should incorporate photoionization by external radiation 
into the Gnat & Sternberg (2007) ionization calculations 
and subsequently follow the ion fractions as the gas cools 
from T ^ 10^ K. This is a computationally elaborate task 
and is beyond the scope of this paper. 

6. ASTROPHYSICAL ORIGIN OF THE BLA - Ovi 
ABSORBER 

In contemplating the nature of this BLA - Ovi ab- 
sorber, we are guided by two key insights, namely (1) the 
absorption is produced by gas at T = 1.4 x 10^ K whose 
ionization is likely governed by both collisional non- 

* The radiative cooling time for a low-density plasma is given by 
tcool = 1.5fcr/(rijjA(T)) { ButSerlaiid & Dopita 1993), where k is 
the Boltzmann's constant and A(T) is the temperature dependent 
radiative cooling efficiency of plasma. For gas at T = 10^ K with 
~ 2% solar metallicity, A ~ 10~^^ erg s"'^ cm'' IGnat fc Sternbere] 
[200W) 
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equilibrium and photoionization reactions (2) the absorb- 
ing structure is part of a large-scale filament of galaxies 
with several luminous galaxies within a ^ 1.5 Mpc 
projected distance. In view of these, a few different sce- 
narios are possible for the origin of the BLA - Ovi ab- 
sorption. The information obtained from probing a pen- 
cil beam region of the absorber may not be sufficient 
to discriminate between these different physical scenar- 
ios. Given the uncertainties encountered in the ioniza- 
tion models, we treat the inferences from the models as 
secondary to what the data is directly suggesting. 

Starburst driven outflows with large velocities from 
galaxies can shock-heat and increase the ionization of 
the surroundi ng gas to produce Ovi (ISh uU fc McKee 
1979; Hcckma n et al 
Oppenheimer fc Pavel 
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IIndebeto"^^^~Shul1 f200l 
Observationally such a 
physical scenario is supported by the detections of Ovi 
systems with [0/H] > 0.2 dex in the Galactic halo as well 
as the extended halos of ot her g alaxies within ^ 500 kpc 
jSavagc ct al. 2002; Fox ct al. '2004'; 'Tumlinson_eLalJ 
[2005: Stockc ct_al. ,2Q06; Wakkcr &._SaYage.,2009). How- 
ever, gas associated with chemical feedback from star- 
formation typically have metallicity higher than the ^ 
— 1.7 dex given by the noneq-CI and photoionization hy- 
brid models for the BLA - Ovi absorber. 

A shock-heated origin for Ovi was suggested by 
iStocke et al.l ([20061 ). In their sample of Ovi absorbers 
and galaxies, 43% (16/37) of the absorbers were as- 
sociated with pairs of Lya clouds within Av = 50 — 
200 km s~^ of each other. Both the Ovi and the Lya 
absorptions [iV(II i) > lO^"^'^ cm~^] were also within 
500 kpc of L* and 100 kpc of 0.1 L* galaxies. About 
half of these Ovi in Lya pairs were also found asso- 
ciated with galaxy groups. The velocity differences be- 
tween the Lya clouds were sufficient to create post-shock 
temperatures for Ovi to be produced. The z = 0.1212 
Ovi absorber towards H 1821 -I- 643 is another example 
in which the hydrogen is seen in a cluster of Lya ab- 
sorbers. In this example, the various Lya components 
were within ~ 100 km s~^ of each other ([Tripp et al.l 
|2001|) . The H i component kinematically aligned with 
the Ovi in that absorber was a BLA with 5 85 km s~^ 
(see Table 3). Several galaxies with impact parame- 
ters ranging from 144 kpc to 3.1 Mpc were also iden- 
tified within I Aw I ^ 500 km s^^ of the absorber imply- 
ing that the sight line was intercepting a galaxy group 
([Tripp et al.l [r998. 2000b). Collisional ionization was in- 
ferred as the likely mechanism for th e ionization in th e 
Ovi - BLA gas phase in that system ([Tripp et al.ll2001[ ). 
The 2 = 0.01027 Ovi absorber is also associated with a 
pair of Lya absorbers separated by |Au| ^115 km s~^. 
Gas clouds moving at those relative velocities, colliding 
with one another, can shock-heat gas to temperatures 
of T, = 1.38 X lO^(w.JlOO)^ ^ 1.8 X 10^ K (assuming 
7 = 5/3 lDraine fc McKeelll993() . consistent with the tem- 
perature for the BLA - Ovi gas phase. Shocks can there- 
fore be one of the mechanisms for the ionization in the 
absorber. 

In a study of Ovi absorbers and galaxies along select 
quasar fields, Chen fc Mulchaey (2009) find some Ovi 
absorbing galaxies with disturbed morphologies indicat- 
ing past events of mergers or interactions. They argue 



that such galaxies likely have a higher cross-section for 
Ovi absorption due to tidally stripped gas contributing 
to the halo gas content. The BLA - Ovi absorber is also 
in a galaxy excess environment 30 galaxies within 
1.5 Mpc projected distance) where the probability 
of galaxy pair interactions is high. Consequently, it is 
plausible that the absorption is from a shock-heated frag- 
ment of gas tidally stripped by NGC 5987 from one of 
the smaller galaxies in the large-scale filament. The oxy- 
gen abundance in the ISM of low-mass irregular galax- 
ies is found to be low (< O.IZ©, ,Skillman_et al. 1989|; 
iKennicutt fc SkillmanI [2001). The tidally stripped in- 
terstellar gas of such a dwarf satellite would be consis- 
tent with the low- metallicity for the BLA - O VI absorber 
given by the hybrid models. The possible difficulty as- 
sociated with this scenario is that tidally stripped gas is 
likely to produce absorp t ion that is kinematically com- 
plex ([Gibson et all [2000; iWakker et"al[[200l iFox et al.l 
[20Tnh . which is not seen in the BLA and Ovi profiles. 

The BLA - Ovi absorption can also be associated with 
the warm- hot intergalactic gas in the galaxy filament. 
As the dominant phase of baryons at z ~ 0, the WHIM 
is expected to be the primary reservoir of gas feeding 
mass into galaxies. The T > 10^ K suggested by the 
data for the absorber corresponds to the temperature 
predicted for the moderately over-dense coo ler compo- 
nent of the WHIM ([Keres fc Hernauistl[2009t) . The BLA 
- Ovi absorber thus presents a unique opportunity to 
probe the interaction of WHIM tied-in with galaxy ha- 
los. The ~ 2% solar metallicity for the absorber, given 
by the hybrid model, is consist ent with the rang e mea- 
sured for the low-z IGM fe.g.. iDanforth fc Shulll [2005). 
Even if the sight line to Mrk 290 is intercepting a nearby 
galaxy's (e.g., NGC 5987) extended halo, the absorption 
could still be from a shock-heated gas cloud accreted 
from the nearby WHIM. Such infiows of matter from 
WHIM are expect ed in the cold gas accretion models of 
galaxy formation (iKeres et al.l[2005HBrooks et al.l 120091 : 
iKeres fc Hernguistl l200'9ir " The presence of paired Lya 
absorbers, and an origin in a galaxy filament, makes the 
absorber an excelle nt candidate for shock-heated WHIM 
([Stocke et al.l[2006[ ). 

7. BLAS WITH METALS AS A TRACER OF WHIM 

In this section we remark on why BLAs with metals 
are an important class of absorption system. Consid- 
erable observational effort has been invested to detect 
warm - hot collisionally ionized gas in the low-z universe, 
as it is predicted to b e the doni inant reservoir phase of 
ordinary matter (e.g., [Bregma^ 2007). The broad-Lya 
absorbers (by definition thermal b{H i) > 40 km s~^) of- 
fer a metallicity iiidependent means to trace w arm - hot 
gas ([Richter et all [20061: IDanforth et al.ll2010[) . However 
the biggest uncertainty associated with the temperature 
measurement in BLAs is in ascertaining the fraction of 
line broadening that is thermal. Blending of closely sep- 
arated absorption components and turbulence within the 
gas can significantly contribute to the broadening of the 
absorption. The extent of this contribution is not read- 
ily assessable from the broad H i profile itself. Detecting 
metals associated with BLAs helps to resolve this prob- 
lem. From the combined 6-values for the metal line and 
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H I, one can solve for the thermal and non-thermal con- 
tributions to the line width and establish a tight con- 
straint on the temperature of the gas. The temperature 
measurement is essential to determine with a certain level 
of precision the ionization fraction of H i and the baryon 
content in BLAs. This value of iV(H) can be used in 
conjimction with single-phase solution to also get an ap- 
proximate measure of the abundance of metals in the 
gas. The absorption system discussed in this paper is an 
example in which such an analysis was possible. 

Metals have been detected in other BLAs. In table 3 
we provide a partial list of BLAs with associated metal 
line detections with |u(BLA) — ^(metal)! < 10 km s~^. 
The close alignment between the H i and metal absorp- 
tion (usually Ovi) can be evidence for the origin of both 
species in a single gas phase. For each example listed 
in table 3, the temperature derived from the 6(H l) and 
6(0 vi) values is T > 10^ K, suggesting a dominant ther- 
mal broadening component for the line widths. The ion- 
ization in the gas phases producing the BLA - Ovi ab- 
sorption in each of these cases could also possibly be un- 
derstood using the nonequilibrium hybrid scenarios sim- 
ilar to the = 0.01027 system described in this paper. 

7.1. Summary 

We have presented HST/COS high S/N observations 
of the z(Ovi) — 0.01027 broad Lya - Ovi absorber in 
the spectrum of the Seyfert 1 galaxy Mrk 290. The ab- 
sorber was discovered by Wakker & Savage (2009) from 
the identification of the Ovi A 1032 A line in the FUSE 
spectrum, which did not cover Lya. The main conclu- 
sions from our analysis are as follows : 

(1) The COS high S/N G130M spectra show Lya ab- 
sorption associated with the z(Ovi) = 0.01027 system. 
In addition, the combined COS G130M and G160M spec- 
tra cover lines from C ii. Si ii, Fe ii, C ill. Si ill. Si IV, 
C IV and Nv. All these ions are non-detections at the 
3 a significance level. The Ovi AA1031, 1037 A lines are 
detected in the FUSE spectrum of Mrk 290. 

(2) The Lya shows absorption from two distinct com- 
ponents with log [N{R i) cm-2] = 14.35 ± 0.01, &(H i) = 
55 ± 1 km s-\ and log [A^(H i) cm^^] = 14.07 ± 0.01, 
6(H i) = 33 ± 1 km s~^ at velocities of = 2 ± 5 km s~^ 
and V = 116 ± 5 km s~^ in the rest frame of the ab- 
sorber. The V = 2 km s~^ component is by definition a 
BLA (b > 40 km s"^). The BLA absorption profile is 
highly symmetric and Gaussian, with a dominant ther- 
mal broadening component of ftthermai/^totai ~ 0.5. 

(3) For the Ovi, we measure line parameters of 
log [iV(Ovi) cm-2] = 13.80 ± 0.05, 6(0vi) = 
34 ± 5 km s-^ The Ovi A1032 A absorption is kine- 
matically well aligned with the BLA, with a Av = 
v{Ovi) - ^(Lya) = -1 ± 5 km s~^. The different b- 
values for the BLA and O vi \delds a kinetic temperature 
of T = 1.4 X 10^ K in the gas, implying that ~ 50% of 
the line broadening is thermal. 

(4) The Mrk 290 sight line is intercepting a filament of 
galaxies at the location of the absorber. The megaparsec 
scale filamentary structure extends over ~ 20 deg in the 



sky, with several luminous galaxies distributed within ~ 
1.5 h^^ Mpc projected distance and Av — 1000 km s^^ 
systemic velocity of the absorber. Among the galax- 
ies, NGC 5987 is closest in impact parameter and ve- 
locity to the absorber. The galaxy is at a projected 
distance of p = 424 kpc, at a systemic velocity of 
\v\ = 71 km s"'^ and oriented almost edge on with re- 
spect to the absorber. NGC 5987 is an Sb galaxy with 
an isophotal major diameter of D25 = 51.7 kpc and a B- 
band magnitude of B = 12.72 or Mb - 5 log /i = -20.4 
corresponding to a luminosity of ~ 2L* . 

(5) Photoionization models do not provide physically 
realistic solutions for the ion ratios in the absorber. The 
model that best satisfies the line strength ratios requires 
the absorbing region to have a very low density of 

6 X 10~^ cm~^ and a path length of L > 3 Mpc. Such a 
large path length would result in Hubble flow broadening 
that is inconsistent with the observed velocity widths 
of H I, Ovi, and their relatively simple and symmetric 
line profiles. Also the temperature prediction of T ~ 
0.6 X 10^ K made by the best-fit photoionization model 
is signiflcantly smaller than T = 1.4 x 10^ K obtained 
from the different BLA and Ovi fe- values. 

(6) We find that the existing simple GIF models, 
noneq-CI cooling models as well as hybrid models of CIE 
and photoionization do not yield a satisfactory explana- 
tion for the ionization in the absorber. The predictions 
most consistent with the observations are from hybrid 
models of non-equilibrium coUisional ionization and pho- 
toionization. Within the frame work of such a hybrid 
model, a metallicity of —1.7 ± 0.1 dex and ionization 
parameter of —1.4 ± 0.1 dex corresponding to a density 
of njj 4 X 10^"'' cm^'', and a baryonic column density 
of A^(H) ~ 4 X 10^'^ cm^^ flts the observed ionic ratios 
at T = 1.4 X 10^ K. The value of the derived metalhcity 
is highly dependent on the validity of the model. In this 
model, the BLA is tracing coUisionally ionized hydrogen 
whereas much of the Ovi is produced via photoioniza- 
tion. 

(7) The galaxy filament and the detection of low- 
metallicity gas at T = 1.4 x 10"'' K suggest that this ab- 
sorber is tracing a shock-heated gaseous structure, con- 
sistent with a few different scenarios including an over- 
dense region of the WHIM in the galaxy filament, or 
the ionized gaseous halo of a galaxy in the filament (e.g. 
NGC 5987). In the latter scenario, the absorber could 
be a fragment accreted by the galaxy from the WHIM. 

(8) Identifying metals associated with BLAs is an effi- 
cient way to detect plasma of 10^''' < 10^ K associ- 
ated with structures that could account for the remain- 
ing ~ 50% fraction of the missing baryons in the low-z 
universe. More such detections of thermally broadened 
Lya and Ovi absorbers are going to be facilitated by the 
high sensitivity spectroscopy capabilities of the Cosmic 
Origins Spectrographs in the far-UV. 
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TABLE 1 
COS Observations of MRK 290 



HST ID 


Grating 


Central Wavelength 


Wavelength Range 


Exposure Duration 






(A) 


(A) 


(sec) 


LB4Q02010 


G130M 


1291 


1136 - 


1430 


964 


LB4Q02020 


G130M 


1300 


1146 - 


1440 


964 


LB4Q02030 


G130M 


1309 


1156 - 


1449 


964 


LB4Q02040 


G130M 


1318 


1165 - 


1459 


964 


LB4Q02050 


G160M 


1589 


1400 - 


1761 


1200 


LB4Q02060 


G160M 


1600 


1412 - 


1773 


1200 


LB4Q02070 


G160M 


1611 


1424 - 


1784 


1200 


LB4Q02080 


G160M 


1623 


1436 - 


1796 


1200 



Note. — All observations were carried out on 28 October 2009 as part of cycle 17 COS-GTO 
proposal ID 1152 (P.I. James Green). 



12 



TABLE 2 
Line Measurements for the 2:(Ovi) 



: 0.01027 Absorber 





I US t r liiiiGut 


Wr 


V 


log [N (cm-2)] 




[~ 


-V, +v] 


Method 






(mA) 


(km s-i) 


dcx 


(km s^-"-) 


(km s -"-) 




Lya 


COS/G130M 




2 ± 5 


14.35 ± 0.01 


55 ± 1 






Fit 








116 ± 5 


14.07 ± 0.01 


33 i 1 








Lycx. 




OZZ It iU 




14.30 ± 0.01 


D4 It 


f- 


220, 70 


AOD 






zyvS it lU 




14.01 ± 0.01 


/II -i- A 


[70, 250 


Lya 




oUo it iU 




14.48 ± 0.01 


OQ -1- A 
96 it 4 


[-220, 250 


AOD 


Ly7 








< 15.2 




[-220, 250 


3 (T 








v{Lya) ~ 2 COMPONENT 










C)\n 1 fi'iO 
kJ VI ±UoZ 






1 ± 5 


13.80 ± 0.05 


9Q -4- 






Fit 


r\\n ^c\'io 

\J VI ±UoZ 




7/1 -1- 1 Q 

/ '4: It lO 


-4 ± 4 


13.82 ± 0.08 


44 IE D 




-75, 75 


AOD 


Ovi 1032 


FUSE 


49 ± 11 




13.63 ± 0.08 


34 ± 4 






AODt 


JNV 1239 


CUb/G130M 


< 35 




< 12.9 






-75, 75 


3 cr 


Nv 1243 


COS/G130M 


< 33 




< 13.0 






-75, 75 


3 a 


C IV 1548 


COS/G130M 


< 48 




< 12.9 






-75, 75 


3 cr 


C IV 1551 


COS/G130M 


< 43 




< 13.1 






-75, 75 


3 rj 


Si IV 1394 


COS/G130M 


< 42 




< 12.4 






-75, 75 


3 a 


C III 977 


GOb/G130M 


< 114 




< 13.3 






-75, 75 


3 (T 


bi III 1207 


GOb/G130M 


< 42 




< 12.1 






-20, 75 


3 a 


C II 1036 


GOb/G130M 


< 36 




< 13.4 






-75, 75 


3 a 


C II 1335 


GOb/G130M 


< 43 




< 13.1 






-75, 75 


3 cr 


bi II liyU 


GOb/G130M 


< 36 




< 12.5 






-75, 75 


3 0- 


■pp nil 

IPC 11 ± A.'-ttJ 


GOS /Gl 301VT 


< 33 




< 13.5 






-75, 75 


3 a 








v(Lya) ~ 116 COMPONENT 










Ovi 1038 


FUSE 


< 31 




< 13.5 






75, 150 


AOD 


Nv 1239 


GOS/G130M 


< 33 




< 12.4 






75, 150 


3 (T 


C IV 1548 


GOS/G130M 


< 39 




< 12.7 






75, 150 


3 cr 


C III 977 


GOS/G130M 


26 ± 27 


103 ± 7 


< 13.2 






75, 150 


3 a 


C III 977 


GOS/G130M 


57 ± 22 


103 ± 7 


13.11 ± 0.22 






80, 130 


3 a 


Si III 1207 


GOS/G130M 


< 37 




< 12.0 






75, 150 


3 a 


C II 1335 


GOS/G130M 


< 40 




< 13.0 






75, 150 


3 a 



Note. — Column 3 lists the rest-frame equivalent width, column 4 gives the velocity ccntroid of line 
components derived from profile fitting, column 5 lists the column density, column 6 is the Doppler parameter 
and column 7 the integration velocity range to derive the respective Wr, N and b. 

All lines other than Lya and Ovi are non-detections at the 3 a significance level. For the two components in 
Lya, different velocity ranges were chosen for measurements of the metal ions. The Ovi AA1031, 1037 A, Ly7 
and C III A977 A are from FUSE. 

For Lya and Ovi A1032 A, we list the line parameters derived from fitting a Voigt function to the pro- 
file and also using the apparent optical depth (AOD) method of Sembach & Savage (1992). The profile fit 
was performed using the Fitzpatrick Sz Spitzer (1997) routine with the fit model for Lya convolved with 
the COS/G130M line spread function described in Ghavamian ct al. (2009). The velocity errors listed are 
dominated by the velocity calibration errors of ~ ±5 km s~^ for both COS and FUSE. For the FUSE 
Ovi A1032 A line a Gaussian kernel of FWHM — 25 km s^^ was convolved with the Voigt model to fit the 
observed profile. 

* The measurement from IWakker fc Savage! S20M) . 
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TABLE 3 
Broad Lya Absorbers with Metals 



Target 


z 


log [N(H l) cm-2] 


6(H I) 


log 


; [Ar(Ovi) cm-^] 


fe(Ovi) 




T 


Reference 






(dcx) 


(km s-i) 




(dcx) 


(km s-i) 


(km s^-"-) 


(K) 




3C 273 


0.00337 


14.27 ± 0.08 


53 + 3 




13.24 + 0.21 


34 + 6 


2 


1.1 X 10^ 


1 


MRC 2251-178 


0.00755 


13.43 ± 0.05 


64 + 4 




13.87 + 0.26 


35 + 5 


8 


1.9 X 10^ 


1 


MRK 876 


0.00312 


14.27 ± 0.08 


77 + 2 




13.18 + 0.21 


29 + 3 


9 


3.3 X 10^ 


1 


FG 1444+407 


0.22032 


13.65 ± 0.05 


86 + 15 




13.94 + 0.07 


36 + 8 


3 


3.9 X 10 


2 


H 1821+643 


0.12147 


13.78 ± 0.17 


or;+37 




14.02 + 0.07 


50+30 


10 


2.4 X 10^ 


2 


H 1821 + 643 


0.26656 


13.63 ± 0.02 


44 + 2 




13.63 + 0.04 


25 + 3 


6 


0.8 X 10^ 


2 


PG 1259+593 


0.31978 


13.98 ± 0.06 


74 + 9 




13.49 + 0.07 


19 + 4 


10 


3.3 X 10^ 


3 










log 


[Af(C III) cm-2] 


6(C III) 








H 1821 + 643 


0.22638 


13.50 ± 0.01 


51 + 2 




12.56 + 0.05 


28 + 4 


1 


1.2 X 10^ 


4 



Note. — A listing of kno wn BLAs with me tals with absorption aligned within |At;| — 15 km s ^ in velocity. The references arc (1) 
iWakker &: Savagel 1120091 ) f2l lTripp et al.l <|2008I 1 (3) Richtcr ct al. (2004) (4) Narayanan, Savage & Wakker (2010). Column 1 is the sight 
line, column 2 the redshift of the absorber, column 3 and 5 arc the column densities of H I and metal ion, column 4 and 6 arc their respective 
6-values, column 7 is the difference in velocity between the Lya and the mctal-line, column S is the temperature derived using the 6-values 
for H I and the metal ion assuming that the two species live in the same gas phase and the listed fc is a quadrature sum of thermal and 
non-thermal b components. 
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Fig. la. — Continuum normalized spectrum of Mrk 290 showing absorption profiles of Lya, and 
Ovi and the wavelength regions of other prominent lines in the z(Ovi) = 0.01027 rest-frame of the 
absorption syste m. The O VI AA1031. 10 37 A, Ly7 and C ll A1036 A, Ly/3 and Ly7 are lower S/N 
FUSE spectra flWakker fc Savagel l2009h and the rest are H ST /COS G130M and G160M observa- 
tions. Every other metal line, except Ovi, is a non-detection at the 3 a significance level. The line 
measurements are listed in Table 2. Ly/3 and Ovi A1038 A are strongly contaminated by Galac- 
tic C II A1036 A and An A1048 A respectively. The wavelength region corresponding to the Ly7 
does not suggest any significant unresolved saturation in H I. Features which are not part of the 
absorption system are labeled with "x" in the respective panels. 
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Fig. lb. — 



Continuation of Figure la. 
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Fig. 2.— To]) Panel : Shows the two component Voigt profile model superimposed on the Lya ab- 
sorption. The vertical axis is continuum normalized flux and the horizontal axis is velocity in the 
rest-frame of the absorber, where v = km s~^ corresponds to 2;(0vi) = 0.01027. Features which 
are not part of the absorption system are marked using "x" in the respective panels. The tick marks 
above the Lya proflle show the centroids of the two absorbing components at u = 2 ± 5 km s~^ and 
V = 116 ± 5 km s~^ respectively. The model Voigt proflles were convolved with the COS/G130M 
line spread function given by Ghavamian et al. (2009). Middle & Bottom Panels : Show the flt to 
the Ovi A1032 A and the contaminated Ovi A1038 A hues respectively. A Gaussian of FWHM 
= 25 km s~^ corresponding to the resolution of FUSE was used to derive the model flt. The flt 
parameters are given in Table 2. 
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Fig. 3.— Apparent column density {Na{v)) profile plots for Lya and Ovi A1032 A. The Ovi A1032 A 
absorption at f = 1 ± 5 km s~^ is kinematically well-aligned with the BLA component of the H I 
absorption at w = 2 ± 5 km s~^, suggesting a single-phase origin for H I and Ovi. The Lya and 
Ovi are binned to approximately the same resolution of ~ 9 km s~^ velocity bins. The Na{v) for 
Ovi A1032 is contaminated at f > 100 km s~^ by a feature unrelated to this system. 
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Fig. 4a. — Top Panel : Galaxies in the foreground field of Mrk 290 within a systemic velocity of u ~ 1000 km s~^ and projected 
distance of ^ 4 Mpc from the BLA-Ovi absorber. The circles represent galaxies and the stars represent sight-lines towards 
AGNs in the field. The distribution pattern of the galaxies shows that the sight line to Mrk 290 is intercepting a megaparsec 
scale filament of galaxies at the location of the BLA - Ovi absorber. Those AGNs with an intervening absorption system 
identified at the redshift of the galaxy filament are plotted in filled star symbols. The size of the circles are proportional to the 
size of the galaxies given by their D25 major diameter. The color coding of the symbols represent various heliocentric velocity 
ranges. The information on galaxies and AGNs in this region of the sky were retrieved from NED. Bottom Panel : A zoom-in 
view of the ~ 1 kpc region centered around Mrk 290. The galaxy IDs are labeled. Among the extended sources, NGC 5987 
is closest in impact parameter (p = 424 kpc) and velocity (|Ai;| = 73 km s~^) to the BLA - Ovi absorber. 
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Fig. 4b. — A Digital Sky Server (DSS) image of NGC 5987 with the line of sight to 
Mrk 290 marked. Figure not included in the astro-ph version due to large file size. 
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Fig. 5 — Photoionization model predictions of column densities for the BLA - Ovi absorber. The 
A^(H l) = 10^^-^^ cm~^ is the column density for the v = 2 km s~^ BLA component of the H I 
absorption. The ionizing background is from the Haardt & Madau (2001) model for z = 0.01027 
and includes UV photons from quasars and young star forming galaxies. The acceptable range of 
column densities for each ion, based on measurement, are highlighted in the photoionization curves 
using thick lines. Except for Ovi, all other ions are non-detections at the 3 a significance level. The 
column density measurements are therefore upper limits for those ions. The multiple X-axes are 
the ionization parameter (log U), total hydrogen column density (A^(H), cm~^), hydrogen number 
density (njj, cm~^), temperature (T, K), gas pressure {p/K, cm~^ K, and size of the absorbing 
region [L = A^(H)/njj, cm). The models assume solar relative elemental abundances of Asplund 
et al. (2009). The metallicity for which the models are able to recover the measured A^(Ovi) is -1.8 
dex of solar. 
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Fig. 6.— CoUisional ionization equilibrium curves for the high ions in the BLA - Ovi absorber based 
on the models of Gnat & Sternberg (2007). The H I column density in these models are set to the 
measured value of 10^"^-'^^ cm~^ (see table 2). We use the updated solar elemental abundances of 



Asplund et al.l (j2009[ ). Except Ovi, the other high ions are non-detections at > 3 cr significance. 
Their column densities are therefore upper limits. The region of the Ovi CIE curve where the 
measured iV(Ovi) ± a[N{Oyi)] is recovered is drawn in thick. Similarly, the regions of the C iv, Nv 
and Si iv CIE curves that are consistent with the measured upper limits of their column densities 
are drawn using thick lines. The vertical dashed line marks T = 1.4 x 10^ K, the temperature given 
by the combined BLA and Ovi 6- values. The vertical dotted line corresponds to the temperature 
upper limit of T = 1.8 x 10^ K if the entire 6(H l) = 55 km s~^ line width for the BLA is thermal. 
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Fig. 7.— Column density predictions from hybrid models where both collisional and photoioniza- 
tion reactio ns are allowed to con tribute to ion fractions. The models were computed using Cloudy 
[ver.C08.00, Ferland et al.l ( 19981 )] . The temperature in the models is fixed to T = 1.4 x 10^ K derived 
for the BLA - Ovi phase, and the A^(H l) = 10^^-^^ cm~^ measured from profile fit to Lya (see table 
2). The acceptable column density ranges based on measurements for the various ions are draw in 
thick. The observed A^(Ovi) is recovered at log U ~ —0.6 for a metallicity of -1.9 dex, assuming 
solar elemental abundances of Asplund et al. ( 20091 ). 
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Fig. 8.— Nonequilibrium coUisional ionization models of Gnat & Sternberg (2007) for 1/10 solar 
( dash - dot line) and so lar {solid line) metallicities and updated solar relative elemental abundances 
of Asplund et al.l ( 20091 ). The acceptable column density ranges for the various ions are drawn in thick 
on each ion curve. Due to rapid rate of cooling, which is higher for the higher gas metallicities, the 
ion fractions of Ovi (and other high ions) are greater than the predictions of CIE in the temperature 
range 10^'^ < T < 10^'^ K. The vertical dashedline marks T = 1.4 x 10^ K, the temperature given by 
the combined BLA and O VI 6- values. The vertical dotted line corresponds to the temperature upper 
limit of T = 1.8 X 10^ K if the entire 6(H l) = 55 km s~^ line width for the BLA is thermal. 
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Fig. 9.— Toy Panel : CIE - photoionization hybrid models computed using Cloudy [ver.C08.00, 
Ferland et al.l (119981 )] for solar metallicity. The Y-axis shows the ion ratio predictions in the hybrid 
models normalized to their corresponding value for CIE. Bottom Panel : Column density predictions 
made by noneq-CI models of Gnat & Sternberg(2007) for a radiatively cooling gas at different tem- 
peratures and a metallicity of 1/100 solar. The vertical dashed line marks T = 1.4 x 10^ K, the 
temperature given by the combined BLA and Ovi 6- values. The vertical dotted line corresponds to 
the temperature upper limit of T = 1.8 x 10^ K if the entire 6(H l) = 55 km s~^ line width for the 
BLA is thermal. To approximately allow for photoionization, the model predictions in the bottom 
panel at T = 1.4 x 10^ K need to be corrected by the values for each ion given in the left panel for 
the appropriate value of log U . 



